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ABSTRACT 

The properties of surface waves in a partially ionized, compressible magnetized plasma 
slab are investigated in this work. The waves are affected by the non-ideal magneto- 
hydrodynamic (MHD) effects which causes finite drift of the magnetic field in the 
medium. When the magnetic field drift is ignored, the characteristics of the wave 
propagation in a partially ionized plasma fluid is similar to the fully ionized ideal 
MHD except now the propagation properties depend on the fractional ionization as 
well as on the compressibility of the medium. 

The phase velocity of the sausage and kink waves increases marginally (by a few 
percent) due to the compressibility of the medium in both ideal as well as Halldiffusion- 
dominated regimes. However, unlike ideal regime, only waves below certain cut-off 
frequency can propagate in the medium in Hall dominated regime. This cut-off for 
a thin slab has a weak dependence on the plasma beta whereas for thick slab no 
such dependence exists. More importantly, since the cut-off is introduced by the Hall 
diffusion, the fractional ionization of the medium is more important than the plasma 
compressibility in determining such a cutoff. Therefore, for both compressible as well 
incompressible medium, the surface modes of shorter wavelength are permitted with 
increasing ionization in the medium. We discuss the relevance of these results in the 
context of solar photosphere-chromosphere. 
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1 INTRODUCTION 

The matter in the Universe is in a partially ionized plasma 
state with varying degree of ionization determined locally 
by ambient thermodynamic conditions. The star-forming 
molecular clouds, accretion discs, solar atmosphere, plan¬ 
etary rings, cometary tail, Earth's ionosphere, and labora¬ 
tory plasma devices are some of the examples of a partially 
ionized medium. The presence of a magnetic field and col¬ 
lision in addition, introduces bewildering variety of scales 
over which dynamics of such a gas manifests itself: as radia¬ 
tion from infalling matter in black holes; as protostellar and 
protoplanetary nebula ensconced in dusty discs; as million 
degree solar corona; as majestic northern lights and finally, 
our own effort to tap the endless source of fusion energy. 
The collision and the magnetic field often have opposite ef¬ 
fect on the plasma dynamics: whereas the magnetic field 
freezes / confines the charged particle motion, collision of¬ 
ten causes diffusion. The neutral gas in such a partially ion¬ 
ized mixture can be acted upon by the magnetic torque via 
collisional momentum exchange with the plasma particles. 
All in all, interplay between the collision and magnetic field 
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may cause relative drift between various constituents of a 
partially ionized gas. To quantify, the relative drift between 
the plasma and neutrals causes ambipolar diffusion while 
relative drift between the electrons and ions causes Hall dif¬ 
fusion. The drift between the electrons and neutrals causes 
Ohm diffusion. Clearly, dynamics of a partially ionized gas, 
in general, cannot be described in the framework of ideal 
magnetohydrodynamics (MHD) and non-ideal effects such 
as Ohm, ambipolar and Hall becomes important. 

In a weakly ionized medium, neglecting plasma inertia, 
a linear relationship between the electric field E and plasma 
current J can be easily derived E = rj ■ J , where rj is the 
diffusivity tensor. However, this approach is not very fruit¬ 
ful in situations where plasma inertia and non-ideal MHD 
effects are simultaneously important i.e. when the plasma is 
not weakly but partially ionized. Typical example includes 
a solar photosphere-chromosphere transition region, Earth's 
ionosphere, protoplanetary discs, discs around cataclysmic 
variables, etc. The collisional, non-ideal MHD effects can 
be included in the dynamics either in the multi-fluid or the 
single-fluid framework. Whereas multi-fluid framework is 
well suited to describe the high frequency fluctuations in 
the medium, the single-fluid framework is often used to de¬ 
scribe the low frequency dynamics of the medium. Given the 
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Figure l. We plot to / Vi n against fractional ionization in 
the above figure. 


complexity of multi-fluid dynamics, it is more than a passing 
curiosity to reduce the multi fluid description to simplified, 
single-fluid, MHD like description. However, the validity of 
single-fluid MHD like description is tied to the fractional 
ionization and ion-neutral collision frequency since the dy¬ 
namical frequency of interest uj must satisfy (Pandey & War- 
dle 2006, 2008) 


( 


UJ 

l^in 


< x 


- 1/2 _ 

e — 



(i) 


where v>i„ is the ion-neutral collision frequency and n e , n n 
are the electron and neutral number densities. The multi¬ 
fluid description (Vranjes et al. 2008; Zaqarashvili et al. 
2012; Soler at al. 2013; Khomenko et al. 2014) is the only 
way out when the frequencies are higher than that given by 
Eq. (1). 

In Fig. (1) tjo/vin is plotted against fractional ionization. 
In Earth's E-region, where X e < 1CV 12 , ultra low frequency 
events are amenable to single-fluid description whereas in 
the solar photosphere-chromosphere where X e < 10 -4 , rela¬ 
tively higher frequency fluctuations can be described in the 
single-fluid framework. Therefore, the quantification high or, 
low frequency depends on the local physical conditions. 

The transition from a weakly to fully ionized region 
poses considerable challenge to the modelling of a partially 
ionized gas. Although such a framework was formally devel¬ 
oped more than half a century ago (Cowling 1957; Braginskii 
1965), the validity of multi fluid framework, Eq. (1) and as¬ 
sociated spatial and temporal non-ideal MHD scales were 
clearly elucidated only recently (Pandey & Wardle 2006, 
2008). For example Hall diffusion becomes dynamically im¬ 
portant at a much lower frequency than the limit derived 
from the fully ionized Hall MHD. Defining pi = m; }Tl m 
as the ion (neutral) mass density where rm >n is the ion (neu¬ 
tral) mass and m is the ion (neutral) number density, 
p = pi + p n as the bulk mass density and uj c i = e B/mi c as 
the ion-cyclotron frequency where e, B and c are the elec¬ 
tron charge, magnetic field and speed of light, respectively, 
the Hall frequency ujh (Pandey & Wardle 2008) 


UJH = — Wei , (2) 

P 

suggests that in a weakly ionized medium (pi —> 0), Hall 
frequency becomes negligible. Therefore, unlike ideal MHD, 
where Hall effect is important only when the dynamical fre¬ 
quency uj is of the order of or, larger than the ion-cyclotron 
frequency, i.e., uj > uJd, in a weakly ionized medium Hall 
effect becomes dynamically important for all frequencies of 
the order of or, above Hall, i.e. uj > ujh which is easily sat¬ 
isfied since ujh ~ 0. Clearly, very low frequency fluctuations 
will be affected by the Hall diffusion of the magnetic field in 
a weakly ionized medium. Therefore, it is not surprising that 
Hall plays important role in the angular momentum trans¬ 
port in protoplanetary discs (Wardle 1999), in destabilizing 
the solar flux tubes (Pandey & Wardle 2012, 2013), in the 
Earth's ionosphere, etc. 

Defining va = B/y/ A-k p as the Alfven velocity in the 
bulk medium, we see that the Hall scale Lh (Pandey & 
Wardle 2008) 


Lh 


( pA 1/2 4 

\p ) ujh ’ 


( 3 ) 


becomes very large in the ujh —> 0 limit, i.e. in a weakly 
ionized plasma Hall operates over a large scale. This feature 
of a weakly ionized medium becomes all the more remark¬ 
able if we recall that in a fully ionized medium Hall scale is 
generally very small (~ ion skin depth) and thus, the Hall 
effect is inconsequential over large scales. 

The MHD waves are capable of carrying energy mil¬ 
lions of miles away from their source of origin and thus, 
play an important role in the space and laboratory plas¬ 
mas. For example, these waves may provide efficient heat¬ 
ing to both fusion (Tataronis & Grossman 1973; Chen & 
Hasegawa 1974; Kapparaff & Tataronis 1977; Mett & Tay¬ 
lor 1992; Pandey et al. 1995) as well as solar coronal plasmas 
(Ionson 1978; Goossens 1994; Parhi et al. 1997b,a, 1998; As- 
chwanden 2009; Goossens et al. 2011). The Alfven waves are 
possible source of turbulence in molecular clouds (Pandey 
& Vladimirov 2007). The past study on wave propagation 
and resonance heating is largely confined to the fully ionized 
medium. Only recently the investigation of wave propaga¬ 
tion in a partially ionized medium has picked up momentum 
due to its importance to space weather (Goodman 1998; 
Khodachenko & Zaitsev 2002; Kazeminezhad & Goodman 
2006; Khodachenko et al. 2006; Pandey et al. 2008; Soler 
at al. 2009, 2010; Goossens et al. 2011; Khomenko & Col¬ 
lados 2012; Zaqarashvili et al. 2012; Pandey 2013; Soler at 
al. 2013; Khomenko et al. 2014). The drift of magnetic field 
through the matter in a partially ionized gas provides new 
pathways through which energy can be channelled to waves 
(Pandey & Wardle 2012, 2013). 

The most important development in the study of MHD 
waves over the last decade in space plasmas has been numer¬ 
ous observational claims of MHD waves in the solar atmo¬ 
sphere [Goossens et al. (2013) and references therein]. This 
has triggered important theoretical development towards ex¬ 
plaining various observations. However, basic MHD model, 
which is employed to explain these observations assume the 
presence of fully ionized matter, a far cry from the state 
of matter below the chromosphere-corona transition region. 
However, since the presence of a magnetic field and den¬ 
sity inliomogeneity considerably complicates the investiga- 
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tion of wave propagation, the neglect of neutral-plasma in¬ 
teraction, although not ideal, is a good starting point. There¬ 
fore on the theoretical front challenge is how to generalize 
the existing MHD wave propagation model to incorporate 
non-ideal, neutral-plasma collision dominated effects. This 
development is crucial for our understanding of underlying 
physical processes at the transition region and thus for the 
resonant heating of the corona and solar wind launching. 

The Sun holds a special place in the investigation of 
MHD waves. The observational validation of various wave 
modes in our backyard allows us to spread our wings to 
distant stars and interstellar medium. Thus it is not at 
all surprising that the large effort has been dedicated to 
understand various wavelike features in the solar atmo¬ 
sphere in the framework of both ideal MHD and Hall MHD 
(Parker 1972; Roberts & Webb 1979; Wentzel 1979; Ed¬ 
win & Roberts 1982; Cally 1985, 1986; Goossens 1994; 
Zhelyazkov et al. 1996; Goossens et al. 2009; Zhelyazkov 
2009). This work will primarily focus on understanding the 
nature of wave propagation in a partially ionized compress¬ 
ible medium and choose as a concrete example solar atmo¬ 
sphere as possible application of the results. Recent investi¬ 
gation of the surface wave in cylindrical filaments (Soler at 
al. 2009, 2010) suggests the important role of non-ideal MHD 
effects for short wavelength fluctuations. We build upon the 
past studies of the surface waves in ideal and Hall MHD and 
generalize it to a partially ionized, compressible medium. 

The paper is organized in the following fashion. The 
basic framework to investigate the surface waves is given in 
Sec. 2. The general dispersion relation for sausage and kink 
mode is derived and various limiting cases are discussed in 
section 3. In Sec. 4 discussion and a brief summary of the 
results is presented and future direction is indicated. 


2 BASIC MODEL 


We shall assume a partially ionized plasma consisting of elec¬ 
trons, singly charged ions and neutral particles. The single- 
fluid MHD-like description of such a gas is given by following 
set of equations (Pandey & Wardle 2008) 


^ + V.(pu) = 0, 

_VP+^ 


(4) 

(5) 


where v = ( pi Vi + p„ v n )/ p, J = n e (vi — v e ) is the current 
density, B is the magnetic field and P = P e + Pi + P n is the 
total pressure. The induction equation is 


dB 

dt 


= Vx 


(v + v B ) X B 



( 6 ) 


where the magnetic drift velocity ( vb ) is defined as 


VB 


(V X B)xb (V X B) 
Vp - o - Vh 


B B ’ ^ 

with b = B/B, J|| = (J ■ b) b , (V X B) ± = V X B — 
— (b-V X B 


b and r/p = pa + p is Pedersen diffusion. 
The Ohm (p), ambipolar (pa) and Hall (ph) diffusivities are 
given as 


p = 


c 

4 7T(J 


PA = 


VA 


andpH = 


v 2 a 

U>H 


( 8 ) 


Here a = cen; (/3 e + Pi) /B is the parallel conductivity de¬ 
fined in terms of plasma Hall parameter 


Pj 


Ucj 

5 

Vj n 


(9) 


which is a ratio between the plasma-cyclotron (uo c j = 
ejB/rruc with e p = ±e) and plasma-neutral ( Vj n = 
PnPnj/pj) collision frequencies. 

To close the above-mentioned set of equations an adia¬ 
batic equation of state P/p 7 = const, will be assumed. Here 
7 is the ratio of two specific heats. Since the background is 
spatially uniform, the spatial derivative of po and Po vanish. 
Here subscript 0 is used to denote background quantities. 
Thus writing pressure and density as / = fo + Sf where the 
fluctuation Sf is much smaller than fo we see that in the 
adiabatic case 


P_ Po ( 5P _ 6p\ Po 

p^pjV+Po 7 poJ pZ’ 


( 10 ) 


and thus we may write SP = c 2 <5p where the sound speed 
Cs = i/v Po/Po- We shall use this equation of state below 
to eliminate perturbed pressure in terms of perturbed den¬ 
sity. Further, equilibrium quantities will be denoted without 
subscript zero below. 

We shall consider a partially ionized, slab of thickness 
2xo having following piecewise constant density and pres¬ 
sure 


p(x) = 


Pin ? Pin if \x\ < xo; 
Pex ■) P ex if |*| > xo , 


( 11 ) 


threaded by uniform magnetic field B = B z parallel to 
the undisturbed surface of tangential discontinuity. Often a 
simplified model of the solar flux tube is modelled by such 
a plasma slab with piecewise constant density and pressure 
(Edwin & Roberts 1982). However, such a discontinuity fun¬ 
damentally changes the behaviour of the Alfvenic vorticity 
propagation in the medium. Whereas in an infinite uniform 
medium, the Alfvenic vorticity in nonzero in the entire vol¬ 
ume, the density jump confines the vorticity to the surface 
layer x = xo only (Goossens et al. 2012). 

We shall further assume that magnetic flux is frozen in 
the plasma slab, moving with v + vb, he. neglect last term 
in the induction equation, (6). This implies that either the 
field aligned parallel current is zero, i.e. Jy = 0 or, the Ohm 
diffusion is unimportant. The linearized equations are 

^+pV-du = 0. (12) 

°*L + v(jJjL + bJb\ = (b^)6b 
dt \ P 4 n p J 4 7r p 


dSB 

dt 


(B ■ V) ( Sv + Svb ) — BV • (Sv + Svb) ■ 


(14) 


Fourier analysing the perturbed quantities as 
exp(—iojt + ikz), and defining L = -j-j — k 2 , andojA = 
k 2 Va > the momentum and induction equations can be 
written in the following form respectively. 



(aP_fc2 c 2) d.7* 
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= v A 
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0 

-k 


0 \ f Sv x \ 
u ) \ Sv v ) 
\ ( SB X /B \ 

J V SB V /B J 


(15) 
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and, 


u> — ir/p L 


0 


ir/H L lo + ik r/p 
_| ,• UJL 


= k 


l UJ 

V A 

0 -1 


5B X /B \ 
SB y /B ) 



(16) 


While writing the above-mentioned equations, we have as¬ 
sumed diffusivities tja and t]h to be constant. Since inves¬ 
tigation of surface wave in a plasma slab is the concern of 
this work, such an assumption will not impact the results. 

Defining a = tv a 2 /tv 2 , fi = c 2 /va and 


A=l~-L,q 2 =Ak 2 
a p 


1 - 


F = 


i Atv 
t\h 


1 + /3 (1 — of) 

. k 2 7] P ' 


1 — a + i 


Fi = 


r)p 


tjh k 2 (1 - a) 


1 — a + i 


. k 2 Tjp 


UJ 


1 + 


y 2 h 
v 2 P j J 


(17) 


Eqs. (15) and (16) can be written in a compact form as 

2 \ / J2 


^ 2 )^ + Ei 


f^2- fc2 )^= 0 » 


2 

dx 2 


d 2 

-k 2 A) 8v x — F 5v v = 0. 


(18) 

(19) 


Note that in the incompressible limit, i.e. when (5 —> oo 
Eqs. (18) and (19) reduces to Eqs. (20) and (21) of Pandey 
(2013) except for a typo i in their F\ in equation (22). 

In the absence of non-ideal MHD effect, i.e. setting r/p = 
r/H = 0, Eqs. (18) and (19) reduces to the following equation 

(^ g2 )^ = °’ (20) 


which is Eq. (5) of Edwin & Roberts (1982) except now 
we are dealing with a partially ionized medium. Clearly, low 
frequency surface wave can propagate undamped in the par¬ 
tially ionized medium (Uberoi & Datta 1998). 

In the absence of Hall diffusion, Eqs. (18) and (19) re¬ 
duces to the following uncoupled equations for 5v x and 8v y 
components 


iaV p L 2 

(l ia fvp) L k 2 

fi UJ 

V pu J 


{uv 2 — uj\ + i k 2 r/p uj ) 5v y = 0 . (21) 

In the incompressible (/3 —» oo) limit, the above-mentioned 
equations are identical to Eq. (55) of Pandey (2013). In 
/3 —>• 0 limit, from Sv x equation above we may write 5v x = 
Qx + Const.. Since fluctuation decays over Pedersen diffu¬ 
sion time-scale (~ l/k 2 r]p), the physical solution requires 
that const. = —Qx o at the surface boundary. Thus, we may 
conclude like an incompressible case that in a purely Ped¬ 
ersen regime the waves will disappear at the surface bound¬ 
ary over ~ 1/k 2 r/p. There is no effect of compressibility on 
5v y equation which also describes wave damping at a rate 
k 2 rip. To summarize, notwithstanding the compressibility 
of the fluid, surface waves are always damped by Pedersen 
diffusion. 

In a homogeneous plasma, we may Fourier analyse the 
a:—dependence as exp(ma:) and the above-mentioned equa¬ 


tions reduce to the following dispersion relation 


+ 1 — a 1 — a + i 


k 2 rjp 


k 2 + a P 


1 — a + i 


+ i 
. k 2 r/p 


k 2 T/p 


1 + 


»7h 

V 2 p 


k z + rr 
n 2 + q 2 / 

= 0 . ( 22 ) 


The dispersion relation can be solved numerically. However, 
we shall analyse it analytically in various limiting cases. 

In the incompressible limit when fi —> oo, we get from 
Eq. (17) A = 1, q 2 = k 2 and the above-mentioned disper¬ 
sion relation reduces to Eq. (23) of Pandey (2013). Since the 
dispersion relation in this limit has been already analysed 
previously, we shall explore fi —> 0 limit i.e. when c 2 -C v\- 
In this limit Eq. (22) becomes 


f k 2 +n 2 \ 
\n 2 + q 2 J 


a 


— a + i 



1 — a + i 


k 2 qp 
UJ 



k 2 qp 
UJ 

No, 


(23) 


which in the absence of Hall (t/h = 0) gives 

(c 2 + ik 2 r/p lo — k 2 v\) (uj 2 + i X Vp w — X va) = 0 , (24) 

where y 2 = n 2 + k 2 . The above-mentioned dispersion rela¬ 
tion, after setting the first bracket to zero, gives mode 


lo ~ ±k va — i 


. k 2 r/p 


(25) 


which describes the damped Alfvenwave. Setting second 
bracket to zero also gives similar formula except now k is 
replaced by x an d tl le damping rate of the wave is ~ x 2 Vp 
and not ~ k 2 r/p. 

In the absence of Pedersen diffusion (r/p = 0), Eq. (22) 
becomes 



2 2 2 

X_ X V h 

k 2 v\ 



(26) 


which for the short wavelength (ujh <S tv a) fluctuations, 
describes the whistler waves in the high frequency {tv a ^ tv) 
limit 

UJ -[ 1 + ^J k2m ■ (27) 

The above-mentioned expression has been derived by bal¬ 
ancing the first and second term in Eq. (26) after noting 
that the last term in square bracket is the dominant term. 
Similarly, electrostatic (VxSE « 0) ion-cyclotron wave in 
the low frequency tv -C tv a limit, 

tv ^ t vh ■ (28) 

is derived by neglecting the first term in Eq. (26)and retain¬ 
ing only the dominant term in the square bracket. 

In the long wavelength limit, i.e. tv a -C uvh, we recover 
usual Alfvencn 2 = tv a and magnetosonic tv 2 = x 2 v\ waves. 
Clearly, wave propagation in limiting cases {fi —> oo and fi —► 
0) display similar behaviour and thus we should anticipate 
that surface waves will also display similar behaviour. 

Although the method of deriving the dispersion relation 
for the surface wave is identical to Pandey (2013), for clarity 
and completeness we shall repeat those steps here as well. 
We start by seeking the solution of Eqns. (18) and (19) as 

8v x — f [exp (—a x) =F exp {a x)] 
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5v y = ih [exp (—a x ) =F exp (a x)] , (29) 

which leads to the following equation for a 

(a 2 — q 2 ) + ( q2 — ^' 2 ) ( a2 — k 2 A) = 0 , (30) 


from where in the incompressible limit we get a 2 = k 2 
and a 2 = k 2 (l — F/k 2 F \) which are same as Eq. (31) of 
Pandey (2013). In order to keep our analytical development 
tractable, we shall approximate the roots of a as 


1 - 


k 2 F\ 


(31) 


which implies that q 2 ~ Ak 2 in Eq. (30) is valid only if 
0 (1 — a) 1. Thus the validity of our analysis requires 
that u}/k > va 0 / {ft — 1) which sets the limit on the com¬ 
pressibility of the medium. 

The four roots in Eq. (31) corresponds to the pair of 
attenuation coefficient (au n i ,cti n2 ) inside and (a ex i ,a ex 2 ) 
outside the slab. Motivated by the fact that planer or cylin¬ 
drical waveguides can support kink and sausage modes, we 
choose the general solution of 5v x and 8v y as a superposition 
of such waves. More precisely, for the sausage wave, inside 
the slab (|a;| < xo) 


sinh(ai„i x) 

5v x (x) = /1 —--- + f 2 

smh(ai n i xo) 

X ( \ ■ t r* sinh(a in i x) 

5vy{x) = 1 f 1 G;„i —--- + 1 f 2 G in 2 

sinh(ai„i xo) 


sinh(aj n2 x) 
sinh(aj„ 2 * 0 ) 
sinh(a in 2 x) 
sinh(a in 2 * 0 ) 


(32) 

(33) 


where 

C , = -, 7 (q 2 i, 2 - Aj k 2 ) 

H Aj [(1 - ap uj + ik 2 r)A] ' 


(34) 


For kink surface-wave, similar expression for the per¬ 
turbed velocities can be given by replacing sinh by cosh. 
The solution outside plasma layer is 

5v x = si e~ aexl (x ~ xo) + s 2 e _ “ ex2 (x ~ xo) , * > *0 


8v x = /3i e“ exl (x+xo) + 0 2 e 1 ** 2 (x+xo) ,x<-x 0 , (35) 

and 


—i8v v = siGexie “ exl(x xo) -| -s 2 G e 


■.2 e 


*ex 2 (x-xq) 


X > Xo 


-iSvy = Pi G ex 1 e“ exl (x+xo) + fa Ge x2 {x+xo) , 

x < —xo ■ (36) 

The knowledge of 5v x and 5v y allows us to calculate the 
perturbed total pressure 


Spr . aui 



1 + f5 (1 — a) d5v x 
1 — pa dx 


+ 


r/p 
a r\ H 


— CL + 


ik 2 ^ k 2 pH ^ 


d5v y 

dx 


(37) 


The electric field components 5E X and 5E y which is required 
for the boundary conditions, can be derived from the gener¬ 
alized Ohm's law 


c8E — — ( 8v + 5vp) X B , 


(38) 


which yields 


c5E x 

B 

and 


= — — Sv , 


C 8Ey 1 

— 5 - =8v x - - 
B a 


Vp_ 

Vh 


1 — a + i 


k 2 q P 


+ i 


. k 2 T)H 


5v v 


(39) 


(40) 


3 DISPERSION RELATION 


We need four boundary conditions across x = xo in order 
to derive the dispersion relation. The first boundary condi¬ 
tion, the continuity of the total pressure across the boundary 
[Spr] = 0 gives 

Xi„ 1 /1 cxini tank (aini * 0 ) + x in2 f 2 ain 2 tanh (cun 2 * 0 ) 

= — (X ex i si a ex i + X ex2 s 2 a ex2 ) (41) 


where Xj =Yj + i Qj Gj and 


Yj = 0- 

Uj 


(1 - aj) 2 Aj 
1 + /3 (1 -aj) 


(42) 


and 



ajPH 



+ * 



(43) 


Since Hall and Pedersen diffusivities in general are a func¬ 
tion of ambient plasma parameters, they will have different 
values inside and outside the plasma slab. However, in order 
to keep the derivation tractable, while deriving dispersion 
relation we have assumed that the diffusivities are constant 
in the plasma. 

By writing the equations in the integral form (as a con¬ 
servation law) for 5v x we obtain fi+f 2 = si + s 2 . Third 
boundary condition is derived by writing the induction equa¬ 
tion in the conservative form 

f TV-17 = 0, (44) 

where 


8U = ( 8v + Sv B ) B — B (8v + Sv B ) , 


(45) 


which gives [ 5U] = 0. The final, fourth boundary condi¬ 
tion is derived by demanding that the continuity of elec¬ 
tric displacement across the surface (Zhelyazkov et al. 1996) 
[8D X \ — 0 where 8D X = K XX 8E X + K xy 8E y with 
2 2 / \ 

I< xx « , K xy « i — . (46) 

v\ v\ \u;hJ 

By imposing the above-mentioned boundary conditions and 
defining Si = Qin/Qe X , S 2 = Si/ (G ex i — G ex 2 ), we arrive 
at the following dispersion relation 


+S 2 


UJA~ 




- 1 


- 1 


“ 2 “ Ti + i -l — (GinlTinl + Gin2Ti n 2) 
L Tin VH 

r fc 2 qp 

(dex 1 T ex 2 Ot eX 2T ex \) % N 


Vh 


UJh 


Vp 


Vh 


— ( - ) ( 1 + — 2 “ ) [GinlTinl + Gin 2 Ti n 2 + S 2 N] — 0 , (47) 
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where 


the following simple form 


Ti — Ti n l + Ti n 2 , 

tanh 


(48) 


Tinl — S3 CXinl ^ coth. J *^* n ) 5 

T in2 = a in 2 ^ ^ ^ {otin2 Xin) , 

T ex 1 = Gi n 2 — G exl/Si — S 3 (Ginl ~ G exl/Si) , 

T ex 2 = Gi n 2 ~ G e x 2 /Sl — S 3 (Ginl ~ Gex^/Sl) , 

N = a exl G exl Tex 2 — a ex 2 G ex 2 Tex 1 • 


with 

0 _ 1 — C Gin2 Si r< [ Id 

‘“ , 3 “ 777 ! ~ 1 ^ 1 — 


1 — C Ginl Si 


A 2 — A-i 


PJ/H J (Pinjpex l) 


Here 


Ai = 


and 


Pin Q e 
Pex Q i 


Cl - 


( W 


V WAi: 


C2 


A 2 = ( Cl 


with 


( w 

ci = - 

V OJH 


P^Ae 


C 2 


(49) 


• (50) 


(51) 


(52) 


1 + 


Vp 

Vh 


+ i 


. Vp 


VH 


C2 = l + t^ 
Vh 


OJ 

(PH 

IP 

U>H 


(53) 


We shall note that when non-ideal MHD effects are com¬ 
pletely absent, i.e. setting vp = Vh = 0 the dispersion rela¬ 
tion Eq. (47) reduces to the following simple form 

(w 2 - u>Ain) +R\ Pt l - {yj 2 - io\ex) tanh(\/)4i K) = 0 , (54) 

V 4±e 

which for Aj = 1 is Eq. (11) of Edwin & Roberts (1982) ex¬ 
cept now it pertains to the partially ionized medium. How¬ 
ever, in the present case due to the compressibility effect, 
Aj V 1 and the above-mentioned dispersion relation is fairly 
complicated. 

There is no general prescription available to solve the 
dispersion relation, Eq. (47). Therefore, we shall analyse it in 
various limiting cases. We first explore the role of Pedersen 
diffusion on the surface waves. To that end, we examine 
the dispersion Eq. (47) in the long wavelength (kx o —> 0) 
limit. Defining R = pex/pin , kx o = K, H = Lp/xo, and 
normalized phase speed Vp = (p/(PAin = we may 

write 


A r . R q Vp 


Qin « ^ Vl , Q e 


,F/Fi « — , 
p Vp 

RQin,Ctt iT-RVp, 


Vh Vh 

Si « 1/R , S 2 ~ — 1 / Gin 2 , S 3 « 1 — /3 , 

Gin 2 « Gex 2 » -t - tanh (« K ) « ^7=^ , 

VpVp VP 

T ex 1 ~ Gin 2 ,7ea;2 « (1 — /3) Gi n 2 ■ 


(55) 


where R q = 1, R for q = in, ex. With the above-mentioned 
approximations, the dispersion relation Eq. (47) reduces to 


VU — (1 — 0) K 2 + i 


(K — 2 i) (3 

HV P 


. y/R-pK (1 — 0) 13 K 


Vp 


2 flV p 2 


(1 + i R) 




(56) 


Note that since in the first bracket, second and third terms 
are the dominant terms, we may balance it with the last 
term in the curly bracket above to yield 


V P « -i 


1/2 


RH 

7 F 


09 - 1), 


(57) 


1 / 2 

which implies that the waves are damping at a rate kVp ■ 
The damping rate is very similar to the well known vis¬ 
cous damping of surface waves in a compressible medium 
(Ruderman 1986, 1991; Ruderman et ah 2000). Clearly, the 
viscosity and Pedersen diffusion plays similar role at the in¬ 
terface. However, detailed physical mechanism in two cases 
are somewhat different. Whereas in the present case it is the 
ion magnetization (determined by a competition between 
the ion-cyclotron frequency against the ion-neutral collision) 
that is responsible for the damping of the wave, it could be 
anisotropic proton viscosity that may be responsible for the 
viscous damping of the waves (Ruderman 1991). 

We shall now analyse the dispersion relation for purely 
Hall case by setting vp = 0. Note that 


— ) - fc 2 - 
A k ~ 


and since 
k 2 




1 - 


£r Fi 
a (3 F 


El = 1 p (PL\ a /P~ a ) ^ x 

a/3F \loh) P(l -a) ' 


in the long wavelength limit, we may write 
Vh F 1 


G,2~ 


Alp Fi (1 — aj) ’ 


(58) 


(59) 


(60) 


for purely Hall case. Defining e = lp/lph and A = 1 + R we 
may write 


C = 


r 2 t /2 

— C ; G 2 = 


As 


S 3 = 


Vp 

Ai 


1 - 


A (yg - i) l 

£ 2 /3 


(61) 


The dispersion relation Eq. (47) becomes 


fAi 


(V£ - 1) + \/ -J- ( R Vp - 1) tanh (sfAi K) 


4 


4-i + 


Aeij-^t^HVA.K) 


= 0. (62) 
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Figure 2. The normalized phase speed of kink wave Vp 
against k xq for ft = 5 (solid line) and /3 = 10 (dashed 
line) is shown in this figure with (frame b) and without 
(frame a) Hall diffusion. 


In the incompressible (P —t oo) limit Eq. (62) becomes 
(yj 2 - u}\ in ) + R (a 2 - wL) tanh(A') 


i0 

uh J 


V 


2 

^ Av 


1 + — tanh(X) 
R 


(1 + R) = 0, 


(63) 


which is Eq. (59) (Pandey 2013) except for a typo in their 
last bracket. 

We shall analyse the dispersion relation, Eq. (62) by 
first noting that since \JWj needs to be positive in order 
for kink wave to propagate in the medium, this implies 
Vp/P < 1. However, since the dispersion relation has been 
derived by assuming /3 (1 — 1/Vf)) 1, combining these 

two conditions gives /3 3> 1. Thus approximating 

= (64) 


and neglecting Vp/ p terms, for a thin slab, the dispersion 
relation becomes 


Vf 


2 P (1 + K) 


R {l + 2PK) +2P + K - 2P H 2 K*/R {1 + R) 


.(65) 


In the absence of Hall {H = 0), the above-mentioned dis¬ 
persion relation reduces to 


Vp ~ 1 + (1 — R) , 


( 66 ) 


which is Eq. (13) of Edwin & Roberts (1982) except here it 
describes the propagation of kink wave in partially ionized 
thin slab. 

We solve Eq. (65) by taking R = 0.25. The choice of 
R is constrained by the equilibrium pressure balance across 
the slab. From Fig. (2) we see that with increasing plasma /?, 
the phase speed of the kink wave also increases. Further, in 
the absence of Hall [Fig. (2(a)] the phase speed propagates 
at all k xo, which is related to the generic nature of the ideal 
MHD like behaviour of the partially ionized medium at low 
frequency (Pandey 2013). However, in the presence of Hall 
[Fig. (2(b)], which incidentally introduces a scale Lh in the 
system, the waves may propagate in the medium only if 


Figure 3. The normalized phase speed of kink wave Vp 
against kx o for Lp = (0.1,0.3,0.5,1) 2+ is shown in 
this figure. 


H 2 K 2 
R (1 + R) 


< R 



1 1 
+ K + 2P ' 


(67) 


It is clear that the cut-off introduced by Hall scale is only 
weakly dependent on the plasma P since Eq. (65) has been 
derived in p 1 limit. 

In thick slab case , when \JAi n K 1, after approxi¬ 
mating \/~A as in Eq. (64), the dispersion relation Eq. (62) 
gives 


Vp 


1 

1 + R — 


H* k 2 ’ 
R 


( 68 ) 


which to the leading order, has no dependence on compress¬ 
ibility. Further, we note that this is Eq. (63) of Pandey 
(2013) except for the absence of 1/R in their Hall term. 
This is due to the fact that the Hall term in their dispersion 
relation Eq. (59) contains 1/(1 + 1/R) rather than 1/(1 +A). 
In the absence of Hall, above equation reduces to Eq. (14) 
of Edwin & Roberts (1982). 

In Fig. (3), we plot Eq. (68) for R = 0.25. It is clear 
that like thin slab, Hall diffusion introduces a cut-off in the 
thick slab as well. Only waves with wavelength larger than 


A > 


2-k L h 

VR (i + R) ’ 


(69) 


can propagate in a thick slab. Clearly, Hall scale, which is an 
artefact of neutral-ion collision in a weakly ionized medium 
(Pandey & Wardle 2006, 2008) introduces such a cut-off. 


4 DISCUSSION AND SUMMARY 


How compressibility does affect the wave propagation in the 
solar photosphere-chromosphere region? To answer this, we 
plot in Fig. (4) plasma P along with the magnetic diffusivi- 
ties by assuming a flux profile 


B = B 0 



(70) 


where Bo = 1200 G. The temperature and density is taken 
from Model C (table 12) of Vernazza et al. (1981). The Ped- 


© 2014 RAS, MNRAS 000, 1 -?? 




















8 B. P. Pandey and C. B. Dwivedi 



Height [Mm] 


Figure 4. The height variation of the Pedersen and Hall 
diffusivities along with the plasma f3 is shown in this 
figure. 

ersen and Hall diffusion profiles have been smoothed by fit¬ 
ting a polynomial to actual model data. It is clear from 
Fig. (3) that in the large part of photosphere, plasma /? 1 

whereas in the upper photosphere (~ 0.5 Mm) and beyond, 
/3 < 1. Clearly, in the middle and upper chromosphere, in¬ 
compressibility is a poor assumption. However, this is only 
valid if the field is 1 kG or, more at the footpoint. For a 
weaker field, e.g. when Bq = 120 G at the footpoint, since 
the value of f5 shown in the figure will be multiplied by 
a factor 100, the plasma can be treated as incompressible 
throughout the entire photosphere-chromosphere. Thus we 
infer from Fig. (2) that the kink waves in thin slab will have 
higher phase speed in the weak field regions. Therefore, in¬ 
tergranular boundaries of the photosphere (known for strong 
vertical fields; (Simon & Leighton 1964)) will launch less en¬ 
ergetic waves than the internetwork regions which often are 
sites of weaker field. It is pertinent to note here that isolated 
intense field have also been detected in the internetwork re¬ 
gions (de Wijn et al. 2009; Sanchez et al. 2010). Thus it is 
quite possible that kink wave of varying energy is excited in 
both network and internetwork regions. 

Note that the Hall diffusion introduces a cut-off, which 
is similar in both thin and thick slabs. Since the Hall scale 
is typically a few km in the photosphere and lower and mid¬ 
dle chromosphere (Pandey 2013), this suggests that waves 
with wavelength larger than 30 — 40 km will propagate in 
plasma slab. The Hall cut-off is not very sensitive to the 
plasma beta. However, Fig. (4) suggests that low-frequency 
sausage and kink wave may not survive in the upper chromo¬ 
sphere altogether. The Pedersen diffusion (which dominates 
Hall in the middle and upper chromosphere) will damp these 
waves whose severity will depend on the wave amplitude and 
ambient plasma parameters. 

We shall note that the present model suffers from 
similar limitations as our previous incompressible model 
(Pandey 2013). Namely, piecewise constant density and pres¬ 
sure profile is an idealization and do not capture the reso¬ 
nant behaviour of the Alfven surface waves (Goossens et al. 
2013). Even with this limitation, the model becomes quite 
complex to handle analytically. Owing to this complexity we 


have analysed only very simple /3 1 case. As noted above, 

plasma beta may become small in the network regions and 
only recourse to investigate the effect of compressibility on 
wave propagation in these regions may be numerical. 

To summarize, wave propagation in partially ionized 
plasma slab only weakly depends on the plasma compress¬ 
ibility. For example waves launched at intergranular bound¬ 
aries may have somewhat larger phase speed than waves 
emanating from internetwork region. However, more than 
the compressibility of the medium, it is the fractional ion¬ 
ization and associated magnetic diffusivities that determine 
the characteristics of normal modes in slabs. 

Following is the summary of this work. 

1. The partially ionized solar photosphere- 
chromosphere plasma is incompressible in the inter¬ 
network region whereas, at intergranular boundaries plasma 
compressibility varies with the scale height. 

2. The surface waves in the compressible plasma 
medium only weakly depends on the plasma (3. The phase 
velocity of sausage and kink waves marginally (by a few 
percent) increases due to the compressibility of a partially 
ionized medium in both ideal as well as Hall diffusion dom¬ 
inated regimes. 

3. Unlike ideal case, in the Hall-dominated regime only 
waves below certain cut-off frequency can propagate in the 
medium. This cut-off for a thin slab has a weak dependence 
on the plasma compressibility whereas for thick slab no such 
dependence exists. More importantly, since the cut-off is in¬ 
troduced by the Hall diffusion, the fractional ionization of 
the medium is more important than plasma compressibility 
in determining such a cutoff. 

4. The ambipolar diffusion always damps long wave¬ 
length fluctuations and is independent of plasma compress¬ 
ibility. 
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